GENERATING SERIES FOR IRREDUCIBLE 
POLYNOMIALS OVER FINITE FIELDS 
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Abstract. We count the number of irreducible polynomials in 
^ . several variables of a given degree over a finite field. The results 

are expressed in terms of a generating series, an exact formula and 
\Q . an asymptotic approximation. We also consider the case of the 

multi-degree and the case of indecomposable polynomials. 
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1. Introduction 



A well-known formula that goes back to Gauss estimate the number 
I n of monic irreducible polynomials among all the N n monic polyno- 
mials of degree n in F 9 [x]: 

(N 

> 4 1 

00 ! N n n 

A proof of this fact from the combinatorial point of view has been 
formalized by Ph. Flajolet, X. Gourdon and D. Panario in [6]. They 
apply this formalism to count irreducible polynomials in one variable 
over ¥ q with the help of generating series. These generating series are 
convergent. The goal of this paper is to count irreducible polynomials 
in several variables. While the formalism is the same, the series are now 
formal power series. At a first glance it seems a major obstruction, but 
despite that the series are non- convergent lots of manipulations are pos- 
sible: computation of the terms of a generating series, approximation 
of this term, Mobius inversion formula,... For comments and examples 
in this vein we refer to the book of Ph. Flajolet and R. Sedgewick [7] 
(and especially Examples 1.19, 1.20, 11.15, Note 11.28 and Appendix A). 

We start by extending the introduction and recall in Section 2 the 
formalism introduced in [6], we consider in Section 3 the case of ir- 
reducible polynomials in several variables of a given degree. Then in 
Section 4 the irreducible polynomials are counted with respect to the 
multi-degree (or vector-degree). Finally in Section 5 we deal with in- 
decomposable polynomials. In each case we give an exact formula, a 
generating series, an asymptotic development at first order and end by 
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an example. References to earlier works are given for each cases after 
the main statements. 

Acknowledgments. I am pleased to thank Philippe Flajolet: this 
article starts after one of his question during a visit at inria Rocquen- 
court, he also explained me how non- convergent series can be nicely 
handled. 



2. Formalism 

Let us recall the formalism of [6]. Let J be a family of primitive 
combinatorial objects and uj 6 X, then = 1 + uj + uj 2 + • • • for- 
mally represents arbitrary sequences composed of uj. The product 
j\f = Y\u>ez it - formally represents all multisets (sets with possible 
repetitions) composed of elements of X. 

In our situation X will be the set of all monic irreducible polynomials 
and Af the set of all monic polynomials. The relation Af = Yluiei ~T^j 
expresses that each monic polynomial admits a unique decomposition 
into monic irreducible polynomials. 

In order to count polynomials over W q we make the following substi- 
tution uj I— > z de&UJ , where z is a new formal variable and degu; is the 
degree of the polynomial uj. We get the generating series: 

that count the total number of monic polynomials and monic irre- 
ducible polynomials. The decomposition of polynomials into irreducible 
factors gives the relation: 

i + n(z) = TT — l - A — = TT - — l — -. 

V ' 11 l _ ^dcgo; 11 Q — z n Y n 

By taking the logarithm we get: 

l0g(l + N(Z)) = V -I n log(l - Z n ) = V 

z — ' z — ' n 

For one variable polynomials we know that the number of monic 
polynomials is N n = q n so that 1 + N(z) = 1 + J2 n >i q n z n = jz^- 
Hence 

1 log; =?, , 

1 — qz n 

by applying Mobius inversion formula to the (convergent) series we get 



it \ i 



n l — qz 
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If we prefer, we can only consider the coefficients of the terms z n in 
Equation (I) to find 

n k n 

k\n k\n 

by applying the classical Mobius inversion formula (/i is Mobius func- 
tion). From this last expression, we deduce that /„ = ^- + O ( 2-^— 1, 
that implies the first statement of the introduction. 



3. Irreducible polynomials: degree 



Let 



bn _ „bn-l 



b„ = I ' y " I and A'„ '' 



v J <? — 1 

Then N n is the number of irreducible, monic, polynomials of degree 
(exactly) n in ¥ q [xi, . . . ,x u ]. The corresponding generating series is 

N(z) = Y J N n z n . 

Let I n be the number irreducible monic polynomials of degree n in 
¥ q [xi, . . . , x u ], we define the corresponding generating series: 

I( Z ) = J2lnZ n . 

Theorem 1. 

(1) ^) = E„^iog(i + iv(^)), 

(2) In = Y,k\n E f[^O g {l + N{z)), 

(3) I n = N n -N 1 - 7V n _! + 0(q b "~^- 2 ). 

In item (3), v and q are fixed while the error term O is for values of 
n that tend to infinity The history for item (3) starts with L. Carlitz 
[3] who proved that almost all polynomials in several variables are 
irreducible. His work has been extended to the case of the bi-degree in 
[4] (as considered in Section 4) and by S. Cohen [5] for more variables. 
Here we get that most of polynomials are irreducible and an estimation 
of the error term. A similar formula has been obtained in [1]: 1 — j^- ~ 

/n — as n grows to infinity. A formula with explicit constants is 



i _ ln_ _ n -n+l 



^ 2q~ n , 



given in [8] for polynomials in two variables 

with n ^ 6. An asymptotic formula of higher order is proved in [11]: 

4 = N n + aiA^-i + a 2 N n _ 2 H h a s N n _ s + 0(iV n _ s _i), where a { are 

constants and the order s is arbitrary large. 

Our proof of (3) appears to be simple compare with these references. 
Finally, in the forthcoming [10], the formula I(z) is also obtained and 
applied to get an approximation: I n = N n —q hn+L '^ 1 (1 + 2/q + 0(l/q 2 )), 
where v and n ^ 5 are fixed, while q grows to infinity. 
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The end of this section is devoted to the proof of Theorem 1. 
Proof. 

Generating series. We follow the formalism of the introduction: 



Hence: 

(2) L(z) : = log(l + N(z)) = £ -I n log(l - z») = £ 

By Mobius inversion formula applied to the series (see comments after 
the proof) we get: 

Another point of view is to make computations term by term: we take 
the coefficients of z n in Equation (2) and get 



[*»]L(*)=x;%- 



k 

k\n 

Then we apply the usual Mobius inversion formula to n[z n ]L(z) and 
find 

(3) I n = J2^[zt]L( Z ). 

k\n 

Asymptotic behaviour. We start with the computation of an ap- 
proximation of N(z) 2 . 

[z n \N{zf = [z n ] [Niz + N 2 z 2 + • • • ) 2 

= 2A n _x • Ni + 0(q bn - 2+b2 - 2 ). 



And for powers k ^ 3 we have: 

[z n ]N{z) k = 0{q b - 2+b2 - 2 ). 
Then we compute the asymptotic behaviour of [z n ]L(z). 
[z n ]L(z) = [z n ] log(l + 

= [z n ]N{z) --[z n ]N(z) 2 + ••• 

= N n - l - (2 • N x • A n _! + 0(q K - 2+b2 - 2 )) + 0(q K - 2+b2 ~ 2 ) 



N n -N 1 -N n _ 1 + 0{q 



b n -2+b2-2\ 
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Now we get: 



k\n 



k\n,k>l 

N n — Ni ■ N n _i + 0(q b - 2+b *- 2 ) + 0(g 6 - 2 ) 
Nn-Ni-N^ + Oiq*"-*-"*- 2 ). 



□ 

Example. For example, in two variables (y — 2), we can compute the 
exact expression of 7 100 . We compute the first 100 coefficients of the 
expansion of L(z) = log(l + N(z)) into powers of z (by derivation an 
expansion of L(z) can be obtained from the inverse of 1 + N(z)). The 
exact value of Jioo is then computed from formula (2) of Theorem 1: 
1 

t 1 ^5151 „5052 „5051 „5050 „4955 , „4953 , r* 4952 , \ 

100 ~ q-V ~ ~ ~ ~ ' 

which has a total of 4385 monomials! When we specialize the ap- 
proximation of Theorem 1 to n — 100 we find that iioo is about 
_L_ (g5i5i _ ^,5052 _ g505i _ ^,5050^ w jth an error term of magnitude g 4954 . 

Finally if we set q = 2 in the exact expression, we get 

Iioo = 4031880625288 . . . 8282220076800, 

a number with 1551 digits. 

Inversion formula. We shall now sketch the proof of inversion for- 
mula for formal power series. Suppose that /, g € with /(0) = 
and 

g{z) = Y,f{z k ) = f{z) + f{z 2 ) + f^) + --- 
k^i 

then 

f(z) = Kn)g{z n ) = ii(l)g(z) + »(2)g(z 2 ) + ^3)g(z 3 ) + • • • 

where /i(n) is the ordinary Mobius function. The proof is in fact the 
same as the classic one: 

EM"M* B ) = EEM")/(^ B )- 

n^l n^l k^l 

In this expression the coefficient of f(z l ) involves only a finite number 
of terms and is X^|jMi) = 1 so that the former sum becomes: 
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Further readings for comments on computations with formal power 
series can be found in [7, Appendix A]. For convergent series, see [13, 
p. 98] and for a general version of inversion theorem, see [14]. 

4. Irreducible polynomials: multi-degree 

Let Nn be the number of irreducible, monic, polynomials of multi- 
degree n— (ni, . . . , n v ) in ¥ q [xi, ...,£„]. By [5] or [11] we have 

q 1 (6 l ,...,5 l ,)e{0,l} v 

The corresponding generating series is 

N(z u ...,z u ) = N ni _ nv -zT---z n u % 

(ni,...,n„)^(0,...,0) 

or for short 

N(z) = J2 N n- 

Let I R be the number of irreducible monic polynomials of degree n 
ia¥ q [xi, . . . , x u ], we define the corresponding generating series: 

n^O 

Theorem 2. 

(1) /U) = E*>i^io g (i + W)), 

(2) 4 = £ fe | gcd(n) ^]log(l + A^ 

(3) ,1,0,. . .,0 ■ JVni,n2-l,n3,...,Ui, ) + 

) if rii^ n 2 ^ ■ ■ • ^ n v . 

For (3) X.-D. Hou and G. Mullen in [11] have already obtained a 
similar formula. The scheme of the proof is the same as in the previous 
section. 

Proof. 

Generating series. The proof is an application of the formalism of 
the introduction, (used in a similar way as in Section 3, see also [4] 
"heuristic proof" of (2.1)). 



i + N(z) = TT - — 



Hence: 

(4) L(z) : = log(l + N(z)) = £ -/„ log(l - = £ 
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where z k = (z* , . . . , z%). Now take the coefficient of z- in Equation (4), 
we get 

k\ gcd(ra) 

By Mobius formula applied to gcd(n)[z-]L(z) we get 

4= E 

k\ gcd(n) 

We can also directly apply Mobius inversion formula to L(z): 

T 

Asymptotic behaviour. We suppose n\ > n 2 ^ • • • > n v . We have: 

fc| gcd(n) 

=4Wa+ E ^b*w*> 

fe| gcd(n),fc>l 



1 , 

=-W„ — ^ (2-^l,0,-,0 ' ^ni-l,n2,...,n v + 2iVo,l,0,...,0 ' ^ni,n2-l,n 3 ,...,n v + ' ' ' 
+ 2iVo,...,o,i • ^m^,... ,^-1,^-1 + O(AT lil)i0i . .. )0 • iV m _i )n2 _i jn3) ... )ni/ )) 
+ O(^2, ,...,0 • ^ ni _ 2 ,n 2 ,...,nJ 

=iV„ — A^i ;0 ,...,0 • ^711-1,712,...,^ + 0(A^o,l,0,...,0 " ^ni,n2-l,n3,..,nj + 
+ O(^2 ) 0,...,0-^n 1 -2,n 2 ,... > nJ. 

Intuitively these formulas express that most of reducible polynomials 
are of type P(xi, . . . , x v ) = {x\ + a) ■ Q(x±, . . . , x v ), where deg Xl Q = 
(deg xi P) — 1, a being a constant. 

□ 

Example. For example, we get J 11)5 = (q 72 — g 67 — g 66 — g 60 + • • • ), 
while iVn j5 = (g 72 — g 66 — g 60 + g 55 ). The approximation of Theo- 
rem 2 specialized to (11, 5) gives that /n i5 is about — ^ (g 72 — g 67 — g 66 ) 
with an error term of magnitude g 61 . Finally if we fix g = 2, we get 
7ii, 5 = 4499945769704095481856 and N u , 5 = 4647462613867219124224. 

5. Indecomposable polynomials 

Definition and known facts. A polynomial P m¥ q [xi, . . . , x u ] is de- 
composable if there exist h e F g [i] with deg /i ^ 2 and Q E ¥ q [xi, . . . , x v ] 
such that 

P(xi, ...,x v ) = ho Q(x 1 , ...,x v ). 
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Otherwise P is said to be indecomposable. For example P(x, y) = 
x 3 y 3 +2xy+l is decomposable (with h(t) = t 3 +2t+l and Q(x, y) = xy). 
Notice that if P — c is irreducible m¥ q [xi, . . . , x u ] for some value cGF, 
then P is indecomposable. 

Indecomposable polynomials have a special interest in regards of 
Stein's theorem and provide families of irreducible polynomials (see 
[15], [12], [2]): 

Theorem (Stein's theorem). If P E F q [xi, . . . ,x„] is indecomposable 
then for all but finitely many values c G ¥ q (the algebraic closure of¥ q ), 
P — c is irreducible in ¥ q [x\, . . . , x v ] . Moreover the number of values c 
such that P — c is reducible is strictly less than degP. 

Let J n be the number of indecomposable polynomials of degree (ex- 
actly) n in F g [xi, . . . , x v ], v ^ 2. Let N n be the number of polynomials 
of degree n (not necessarily monic) in ¥ q [xi, . . . , x u \. 

N n = q K -q K -\ 



We define the following Dirichlet series: 

^— ' n s n s n s 

nJSl n^l n^l 

We also define a kind of Mobius function /i(d,n): 

Kd,n)= £ (-lM+-+^, 

d=di | d,2 1 c?3 ■ • • | d k | d k + x =n , 
di<rf2<---<rffc<™ 

where the sum is over all the iterated divisors of n (k not being fixed). 
We define fi(n,n) = 1. For a fixed n, we define 1 < £ < £' to be the 
first two divisors of n. 

Theorem 3. 

(1) J(s) = §}, 

(2) J n = E d |n,l^n^( d ' 71 )^ 

(3) 7/n zs £/ie product of at least three prime numbers: 

J n = N n - q ^N nj+ 0(q w '- 2 N f ). 

The item (3) has already been obtained in [2] and [9], where some ex- 
plicit constants are computed. For example in two variables, and when 
n is the product of at least three prime numbers, Theorem 5.1 of [2] 

yields the same first order estimate of J n : 

n^—nrlSlkz - n r ■ 
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Proof. 

Generating series and inversion formula. The decomposition of 
a decomposable polynomial is unique once it has been normalized, see 
[2, Section 5]: 

Lemma 4. Let P = h o Q be a decomposition with degh ^ 2, Q 
indecomposable, monic and its constant term equals zero. Then such a 
decomposition is unique. 

By induction it implies the following lemma (see [2]): 

Lemma 5. 

J n = N n - q*- 1 x J d . 

d\n,l^.d<n 

From Lemma 5 we get: 

N n = Yl ^ x Jd - 

d\n,lS^.di^n 

From the point of view of Dirichlet series it yields: 

N(s) = J(s)-F(s). 
Then by induction we also deduce from Lemma 5 that 

d\n,l^.d^.n 

Asymptotic behaviour. From this formula we deduce: 

J n = /i(n, n)N n + /2 (n/£, n) N* + O (jj, (n/£', n) N^j 

= N n -q^N^+0((q W '- 2 -q^)N, 

= N n - ^Nn + O (q W '- 2 N f ) . 



□ 

Example. For example, in two variables, we have iV 100 = g 5151 — g 5050 
and we compute J 100 by formula (2) of Theorem 3. We get J wo = 
g5i5i _ ^,5050 _ ^1327 _|_ ^,1276 _ ^354 _| _ While the approximation given 

in Theorem 3 is that J wo is about g 5151 - g 5050 - q 1327 + g 1276 with an 
error term of magnitude g 355 . 
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